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Figure 1: Diagrams illustrating one- and two-body jEFT EM currents entering at LO (e Q^^), NLO (e Q ' ), 
N2LO (eQ^), and N3LO (e Q ' ). Nucleons, pions, and photons are denoted by solid, dashed, and wavy lines, 
respectively. 



In this contribution, we present a set of Green's function Monte Carlo (GFMC) calculations of 
magnetic moments (m.m.'s) and Ml transitions in A < 9 nuclei, which has been recently reported 
in Ref. In these calculations, nuclear wave functions (w.f.'s) are constructed from a Hamil- 
tonian consisting of the Argonne-vig two-nucleon [||] and Illinois-7 three-nucleon potentials |^], 
with which the computed GFMC ground- and excited-state energies are found to be in good agree- 
ment with experiments The electromagnetic (EM) current operator includes, in addition to the 
standard one-body convection and spin-magnetization terms for individual protons and neutrons, a 
two-body meson-exchange-current (MEC) component. The latter is constructed within two distinct 
frameworks, namely the standard nuclear physics approach (SNPA) illustrated in Refs. [Q, |5|], and 
the pionfull chiral effective field theory (;^EFT) formulation of Refs. [||,^|8|]. In what follows, we 
summarize on the methods and results discussed in Ref. 

1. GFMC Method and the Nuclear Hamiltonian 

The EM transition matrix elements are evaluated in between w.f.'s which are solutions of the 
Schrodinger equation, with a nuclear Hamiltonian, H, consisting of a kinetic term plus two- and 
three-body interaction terms — in the present case, the Argonne-vi8 and Illinois-7, respectively. 
Nuclear w.f.'s are constructed in two steps. First, a trial variational Monte Carlo w.f. (*Pr), which 
accounts for the effect of the nuclear interaction via the inclusion of conelation operators, is gener- 
ated by minimizing the energy expectation value with respect to a number of variational parameters. 
The second step improves on by eliminating excited states contamination. This is accomplished 
by the GFMC calculation which propagates the Schi^odinger equation in imaginary time (t). The 
propagated w.f. ^{t) = for large values of T, converges to the exact w.f. with eigen- 

value Eq. Ideally, the matrix elements should be evaluated in between two propagated w.f.'s. In 
practice, we evaluate mixed estimates in which only one w.f. is propagated, while the remaining 
one is replaced by ^'t- The calculation of diagonal and off-diagonal matrix elements is discussed 
at length in Ref. [Q] and references therein. 
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The nuclear EM current operator is also expressed as an expansion in many-body operators. 
The current utilized in the calculations accounts up to two-body effects, and is written as: 

j(q) = D'(q)+Ij'7(q)> d-D 

where q is the momentum associated with the external EM field. The one-body operator at lead- 
ing order, i.e. the impulse approximation (lA) operator, consists of the convection and the spin- 
magnetization currents associated with an individual nucleon [[T]], and it is diagrammatically repre- 
sented in panel (a) of Fig. [1} 



2. ;t:EFT and SNPA EM currents 

In the calculations, two models for the EM two-body MEC operators are tested, namely the 
pionful ;^EFT and SNPA models. The ;^EFT current operators are expanded in powers of pi- 
ons' and nucleons' momenta, Q, and consist of long- and intermediate-range components which 
are described in terms of one- and two-pion exchange contributions, as well as contact currents 
which encode the short-range physics. These operators involve a number of Low Energy Constants 
(LECs) which are then fixed to the experimental data. Currents from pionful ;^EFT including up 
to two-pion exchange contributions were derived originally by Park, Min, and Rho in covariant 



perturbation theory [10]. More recently, Kolling and collaborators presented EM currents obtained 



within the method of unitary transformations [ |11| , |12| ]. Here, we refer to the EM operators con- 
structed in Ref. [^, ^,|8|], in which time-ordered perturbation theory is implemented to calculate the 
EM transition amplitudes. These EM operators are diagrammatically represented in Fig. |I| where 
they are listed according to their scaling in eQ, (where e is the electric charge). 

Referring to this figure, one-body contributions enter at LO, panel (a), and N2L0, panel (d), 
and they are the lA current operator at LO and its relativistic correction, respectively. The NLO 
term involves seagull and in-flight long-range contributions associated with one-pion exchange 
(OPE). At N3L0 we include the two-pion-range contributions of diagrams (e)-(i), the one-pion- 
range tree-level current involving a ynNN vertex of order eQ^, diagram (j), the contact cuixents of 
diagram (k), as well as the one-loop corrections of diagrams (l)-(o). The two-body operators have 
a power-law behavior at large momenta, therefore a regularization procedure is implemented via 
the introduction of cutoff function of the form exp{-Q'^/A'^) where A = 600 MeV. 

The contact currents of diagram (k) are of minimal and non-minimal nature. The former are 
linked to the xEPT potential at order via current conservation; therefore they involve the same 
LECs entering the ;^EFT NN interaction, and can be taken from fits to the NN scattering data. We 



use the values obtained from the analysis of Refs. []16[], with cutoff A = 600 MeV. Non-minimal 
LECs entering the contact and tree-level currents at N3L0 — diagrams (j) and (k), respectively — 
need to be fixed to EM observables. The fitting procedure has been implemented by Piarulli et al. 
in Ref. [^]. In that work, LECs multiplying isovector operators in the tree-level current are fixed by 
saturating the A-resonance JTO| ] (a common strategy adopted, for example, in Refs. [13, 14, 15]). 



The remaining thi^ee LECs ai^e fixed so as to reproduce the deuteron, ^He, and m.m.'s. 

The second model for the EM MEC operators utilized in the calculations is the SNPA model. 
Two-body currents in the SNPA formalism, described at length in Refs. [^, ^ and references 
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Figure 2: Left: Magnetic moments in nuclear magnetons for A < 9 nuclei. Black stars indicate the ex- 
perimental values JitI [T^ , while blue dots (red diamonds) represent preliminary GFMC calculations which 
include the lA one-body EM current (full jEFT current up to N3LO). Predictions are for nuclei with A > 3. 
Right: Magnetic density in nuclear magnetons per fm^ for selected nuclei, including only the lA current 
contribution. 



therein, are separated into model-independent (MI) and model-dependent (MD) terms. The former 
(MI) ai^e derived from the AW potential, and their longitudinal components satisfy, by construction, 
current consei"vation with it, thus their short-range behavior is consistent with that of the potential. 
The dominant terms, isovector in character, originate from the static part of the potential, which is 
assumed to be due to exchanges of effective pseudoscalar (PS or "TT-like") and vector (PV or "p- 
hke") mesons. The associated currents are then constructed by using the PS and PV propagators, 
projected out of the static potential At large inter-nucleon separations, where the AW potential 
is driven by the OPE mechanism, the MI current coincides with the standard seagull and pion-in- 
flight OPE currents diagrammatically illustrated in panels (b) and (c), respectively, of Fig. [I[ The 
MD currents are purely transverse, and unconstrained by current conservation. The dominant term 
is associated with excitation of intermediate A isobars. Additional and small MD currents arise 
from the isoscalar p;r7and isovector w Try transition mechanisms [||, ||]. 



3. Results 

Results for the m.m.'s of A < 9 nuclei indicate that two-body MEC corrections evaluated in 
both the SNPA and ;tEFT models are qualitatively in agreement, and, when large, they boost the 
LA in the direction of the experimental data. We summarize the m.m.'s calculations in the left panel 
of Fig. ^ where we show the results obtained with the ;^EFT model. In this figure, black stars 
represent the experimental data [17 , 18] — there are no data for the m.m. of ^B.' For completeness, 
we show also the experimental values for the proton and neutron m.m.'s, as well as their sum, 
which corresponds to the m.m. of an S-wave deuteron. The experimental values of the A = 2-3 
m.m.'s have been utilized to fix the LECs, therefore predictions ai^e for A > 3 nuclei. The blue dots 



'Electromagnetic moments of lithium isotopes have been most recently measured in Refs. pof . At the present 
time, the tables and figures of this contribution and the preprint of Ref. ^ show, for these nuclei, the experimental data 



taken from Refs. [0,|T8|l. 
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Figure 3: Left: Ratio to the experimental Ml and E2 transition widths in A < 9 nuclei. Black stars with error 
bars indicate the experimental values [ [l7[ |l8[, while blue dots (red diamonds) represent GFMC calculations 
which include the lA one-body EM current (total jEFT current up to N3LO). Right: Ml transition density 
in nuclear magnetons per fm^ for selected nuclei, including only the lA current contribution. 



labeled as GFMC(IA) represent lA theoretical predictions. The GFMC(IA) results reproduce the 
bulk properties of the m.m.'s of the light nuclei considered here. In particular, we can recognize 
three classes of nuclei, that is nuclei whose m.m.'s are driven by an unpaired valence proton, or 
neutron, or 'deuteron cluster' inside the nucleus. This behavior can be appreciated by looking at 
the lA magnetic densities represented in the right panel of Fig. |2[ where the red upward-pointing 
triangles are the contribution from the proton spin, jUp[Ppt(r) — Ppi{r)], and similarly the blue 
downward-pointing triangles are the contribution from the neutron spin. The green diamonds are 
the proton orbital (convection current) contribution, and the black circles are the sum. For example, 
we can see that the m.m.'s of 'Li and ^Li are driven by the unpaired proton, while the m.m. of ^Li 
it is due to a combined effect of the unpaired neutron acting against the proton. 

In the left panel of Fig. ^, predictions, which include all the ;^EFT EM current contributions 
illustrated in Fig. |l], are represented by the red diamonds labeled GFMC(TOT). In all of the cases 
considered here — except for ^Li and ^Be for which the I A results are already in very good agree- 
ment with the experimental data, the predicted m.m.'s are closer to the experimental data when 
the MFC corrections are added to the lA results. MFC corrections are particularly pronounced 
in the isovector combination of the A =9, T = 3/2 nuclei's m.m.'s, for which the MFC SNPA 
(;fFFT) correction provides ~ 20% (~ 30%) of the total calculated isovector contribution. While 
the SNPA and ;^FFT models are in a reasonable good agreement when predicting the isovector 
m.m.'s — which are driven by the long-range NLO OFF contribution, we find that isoscalar m.m.'s 
evaluated within the ;^FFT model are usually in a better agreement with the experimental data |jl|]. 

In the left panel of Fig. |3[ we show FM ti^ansitions induced by the Ml and F2 operators in 
A <9 nuclei — F2 transitions are provided in lA only. In this figure, we show the ratios to the 
experimental values of the widths [17, The latter are represented with the black stai^s along 
with the associated experimental error bars, while the GFMC(IA) and GFMC(TOT) predictions are 
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again represented by blue dots and red diamonds, respectively. For the Ml transition in lA, we also 
provide their transition densities which are illustrated in the right panel of Fig. ^ As before, the 
red upward-pointing triangles are the contribution from the proton spin term, the blue downward- 
pointing triangles are from the neutron spin, the green diamonds are from the proton orbital tenii, 
and the black circles are the total lA contribution. For example, for the lithium isotopes, the Ml lA 
transitions are predominantly from the proton spin term, i.e., these are almost pure proton spin-flip 
transitions. While, for ^Be and ^B, the neutron spin term is the most important, but with some 
contribution from the proton spin and orbital terms. The Ml results summarized in the left panel 
of Fig. |3[, indicate that, also for these observables, predictions which account for MFC corrections 
are closer to the experimental values, but for the transition in *^Li, for which the experimental error 
is large, we cannot determine whether the GFMC(TOT) prediction is a better one. 

This work is supported by the U.S. Department of Fnergy, Office of Nuclear Physics, under 
contracts No. DF-AC02-06CH1 1357 and No. DF-AC05-06OR23177, under the NUCLFI 
SciDAC-3 grant, and under NSF grant PHY- 1068305. 
References 

[1] S. Pastore, S. C. Pieper, R. Schiavilla, and R. B. Wiringa, arXiv: 1212.3375. 
[2] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev C 51, 38 (1995). 
[3] S. C. Pieper, MP Conf. Proc. 1011, 143 (2008). 

[4] L. E. Marcucci, M. Viviani, R. Schiavilla, A. Kievsky, and S. Rosati, Phys. Rev C 72, 014001 (2005). 

[5] L. E. Marcucci, M. Pervin, S. C. Pieper, R. Schiavilla, and R. B. Wiringa, Phys. Rev C 78, 065501 
(2008). 

[6] S. Pastore, R. Schiavilla, and J.L. Goity, Phys. Rev. C 78, 064002 (2008). 

[7] S. Pastore, L. Girianda, R. Schiavilla, M. Viviani, and R. B. Wiringa, Phys. Rev. C 80, 034004 (2009). 

[8] M. Piarulli, L. Girianda, L.E. Marcucci, S. Pastore, R. SchiaviUa, and M. Viviani, Phys. Rev. C 87, 
014006 (2013). 

[9] M. Pervin, S. C. Pieper, and R. B. Wiringa, Phys. Rev. C 76, 064319 (2007). 
[10] T.-S. Park, D.-P Min, and M. Rho, Nucl. Phys. A596, 515 (1996). 
[11] S. KolUng, E. Epelbaum, H. Krebs, U.-G. Meissner, Phys. Rev. C80, 045502 (2009). 
[12] S. KolUng, E. Epelbaum, H. Krebs, and U.-G. Meissner, Phys. Rev C 84, 054008 (201 1). 
[13] Y.-H. Song, R. Lazauskas, T.-S. Park, and D.-P Min Phys. Lett. B 656, 174 (2007). 
[14] Y.-H. Song, R. Lazauskas, and T.-S. Park, Phys. Rev. C 79, 064002 (2009). 
[15] R. Lazauskas, Y.-H. Song, and T.-S. Park, Phys. Rev. C 83, 034006 (2011). 

[16] D.R. Entem and R. Machleidt, Phys. Rev. C 68, 041001 (2003); R. Machleidt and D.R. Entem, Phys. 
Rep. 503, 1 (2011). 

[17] D. R. Tilley, C. M. Cheves, J. L. Godwin, G. M. Hale, H. M. Hofmann, J. H. Kelley, C. G. Sheu, and 
H. R. Weller, Nucl. Phys. A 708, 3 (2002). 

[18] D. R. Tilley, J. H. Kelley, J. L. Godwin,D. J. Millener, J. E. Purcell, C. G. Sheu, and H. R. Weller, 
Nucl. Phys. A 745, 155 (2004). 

[19] D. Borremans et ai, Phys. Rev. C 72, 044309 (2005). 

[20] R. Neugart et ai, Phys. Rev. Lett. 101, 132502 (2008). 



6 



